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Abstract 

The model of p Ising spins coupled to 2d gravity, in the form of a sum over 
planar (jr' graphs, is studied and in particular the two-point and spin-spin correlation 
functions are considered. We first solve a toy model in which only a partial summation 
over spin configurations is performed and, using a modified geodesic distance, various 
correlation functions are determined. The two-point function has a diverging length 
scale associated with it. The critical exponents are calculated and it is shown that all 
the standard scaling relations apply. Next the full model is studied, in which all spin 
configurations are included. Many of the considerations for the toy model apply for 
the full model, which also has a diverging geometric correlation length associated with 
the transition to a branched polymer phase. Using a transfer function we show that 
the two-point and spin-spin correlation functions decay exponentially with distance. 
Finally, by assuming various scaling relations, we make a prediction for the critical 
exponents at the transition between the magnetized and branched polymer phases in 
the full model. 



1 Introduction 

The work of Knizhnik, Polyakov and Zamolodchikov (KPZ) Q as well as of David, Distler 
and Kawai (DDK) Q made it possible to understand many aspects two-dimensional 
quantum gravity coupled to conformal field theories. At the same time, it became clear 
that the models of dynamical triangulations coupled to matter fields provide us with a 
statistical mechanical realisation of the models described by KPZ and DDK, in the same 
way as the two-dimensional Ising model at its critical point can be viewed as a conformal 
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c = 1/2 field theory. Two aspects of 2d quantum gravity coupled to matter fields remain 
puzzling: the c = 1 barrier and the concept of correlation length. The formulae for critical 
exponents derived by KPZ cease to be valid for c > 1 and the critical exponents are derived 
by general scaling arguments applied to globally defined operators. At no point is the 
concept of a divergent correlation length introduced. For ordinary statistical systems the 
divergence of a correlation length when a critical temperature is approached is believed 
to be the underlying reason that general scaling arguments work well. The difficulties of 
defining a local length scale in quantum gravity are well known and have so far prevented 
a proper treatment of correlation functions by means of continuum methods. 

Working entirely in the context of dynamical triangulations the problems mentioned 
are not seen directly. Statistical models with multiple Ising spins living on dynamically 
triangulated surfaces are perfectly well defined even if c > 1 and they have a critical point. 
At least superficially, the only difference compared to c < 1 is that we cannot solve the 
theory. Nevertheless, low temperature expansions and mean field calculations seem reli- 
able for c — > oo, as is confirmed by the agreement between the theoretical calculations and 
Monte Carlo simulations of the systems with large c |5|. The picture which emerges for 
large c from the low temperature expansion [Q] is as follows: for large /3 (low temperature) 
there is a magnetized phase for which ^str = separated from a branched polymer 

phase, where ^str = 1/2, by a transition at a critical (3* where 7*^^, = 1/3. 

Likewise it has been possible by the use of dynamical triangulations to address the 
question of correlation length in two-dimensional quantum gravity. A two-point correlator 
between "punctures" has been calculated in pure gravity as a function of a "quantum" 
geodesic distance and it is found that standard scaling relations, known from the theory of 
critical phenomena, are satisfied, although with unusual critical exponents So far it 
has not been possible in 2d quantum gravity to calculate correlation functions of matter 
fields as a function of the geodesic distance. The attempts to measure the spin-spin 
correlation functions by Monte Carlo simulations and to define a divergent correlation 
length as one approaches the critical point, have so far been ambiguous |^. The 
question arises as to whether there is a divergent spin-spin correlation length associated 
with the phase transition between a magnetized and a non-magnetized phase in the two- 
dimensional Ising model coupled to gravity. Rather surprisingly from a continuum point 
of view, we can, using dynamical triangulations, begin to answer this question in the limit 
of large c, i.e. in the limit where a large number of Ising models is coupled to quantum 
gravity. 

The statistical model we will define and solve in the following sections is a toy model 
of 2d quantum gravity coupled to Ising spins in the sense that the correct summation 
over all triangulations is performed, but not all the spin configurations are included. The 
spin configurations we include are precisely the ones which dominate in the large c limit, 
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at large /3, namely those for which the domains are connected in a tree-like fashion with 
domain boundaries of minimal size. This model allows us to calculate the two-point 
functions (using a certain definition of distance, which is based on the geodesic distance) 
and extract the critical exponents. We will verify that these exponents satisfy standard 
scaling relations and that the geometrical interpretation of some of the exponents is 
related to the fractal structure of the underlying "space-time" , in agreement with general 
arguments; that is, we explicitly verify that the exponent v is related to the Hausdorff 
dimension hy v = l/dff. The model has a third order transition between a tree- like (i.e. 
branched polymer) phase and a magnetized phase. It is found that there is a diverging 
correlation length in the tree phase associated with the geometric two-point function, but 
no diverging correlation length associated with the spin-spin correlation function. 

Next we turn our attention to the full model of p independent Ising spins coupled 
to 2d gravity, which has a central charge of c = p/2. Following the analysis of the toy 
model, we show that various two-point functions decay exponentially with distance in 
the magnetized phase. Again there is a diverging correlation length associated with the 
geometric two-point function. By assuming that all the standard scaling relations still 
hold for this model and making a few further fairly modest assumptions, we show that 
the critical exponents for the transition between the magnetized and branched polymer 
phases, in the full model, are the same as those in the toy model. 

The rest of this article is organized as follows: in section § we define the toy model, 
define a variant of geodesic distance in the model and calculate the two-point function. 
In sections |3| and ^ the critical exponents and Hausdorff dimension are calculated for the 
toy model, whilst section ^ discusses the spin-spin correlation function. In section ^ we 
verify by direct calculation in an external magnetic field, that the magnetic exponents 
found by scaling arguments are indeed the correct ones. In section |^ we address some of 
the questions mentioned above for the full model, i.e. in the model where we sum over all 
the spin configurations. Finally, section |8| contains our conclusions. 

2 Toy model, without a magnetic field 
2.1 Definition 

Before looking at the full model of p spins coupled to 2d gravity (in the form of a sum over 
planar (f)^ graphs), we will solve a toy model, which is very similar to the one studied in Q. 
As for that model we will show that there is a magnetized phase for which jstr = and 
a branched polymer phase with ^str = ^- On the boundary 7*^^ = i, which is the same 
value as that in Q and agrees with the result 7*^,, = 7/(7 — 1) for 7 = — |, the pure 
gravity value, which one might expect from the calculation in [^. 
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For our model we sum over all possible rooted planar (jr" graphs (let us denote this 
set of graphs by Ql), but we only sum over a subset of the possible spin configurations, 
namely those for which the domains are connected in a tree-like fashion, with at most 
one link connecting any two domains. In this paper we shall only consider planar graphs, 
that is, X = 2 throughout. Each vertex will be weighted with a factor of x and has p 
independent spins on it which can take the values ±1 (i.e. vertex i has spins Sf on it, 
with Q = 1, • • • Links joining vertices with dissimilar spin configurations will be given 
a factor of e~^^ for each spin flavour which differs. Thus the grand canonical partition 
function, T, is 

T= E^^E n«-pfE/3({^f^"-i))' (1) 

Gagi {sy <ij> \a=l / 
where the first summation is over graphs, with N being the number of vertices in the 
graph. The product is over the nearest neighbour pairs on the graph G (referred to 
as "links") and the second summation is over the following set of spin configurations. 
Take the graph G and decompose it into a set of one-particle irreducible ("IPI") graphs 
(which we shall call "blobs") connected in a tree-like fashion (see fig. |l|); note that these 
blobs are essentially just minimum necked baby universes ("minbu"s). The blobs are fully 
magnetized, that is, all the vertices within a IPX blob have the same set of spins, and we 
will sum over the 2^ possible spin configurations for each blob (except the root blob for 
which all the spins will be fixed to be +1). 

Figure 1: Measuring distances 
Distance 1 2 3 4 




In the next sections we will solve this model exactly and also calculate various corre- 
lation functions, but in order to do this we need first to define how distances are to be 
measured on the graphs. One possible measure of the distance between two vertices A 
and B, is the geodesic distance between them. That is, find the shortest path from A 
to B along links, counting a distance of one unit for every link traversed — this gives the 
geodesic distance between the two points. Unfortunately, this is quite difficult to deal 
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with analytically, although some results are known for the pure gravity case |^. In this 
paper we use a slightly different definition of distance, which is much simpler to handle. 
The distance between two vertices will be defined as the shortest path between them along 
links, counting a distance of one unit only for links that separate two IPI blobs. 

Since we are interested in correlations between the spins in the root blob, which are 
held fixed and spins further away, distances will be measured from the root blob. Thus all 
the vertices in the root blob have a distance zero, all the vertices within blobs connected 
by a single link to the root blob are at distance one and so on (see fig. |l|). It will turn out 
that the average number of vertices in each blob is very small, so that for many graphs 
this definition of distance will be quite similar to the geodesic distance, especially in the 
branched polymer phase for which the partition function is dominated by tree-like graphs. 

The advantage of measuring distances in this way is that it makes it very easy to 
define a transfer function, f{y), 

fiy) = Z" (x(l - A2/)"i) + ^y^ (2) 

where A = 2x ^1 + e^^'^^^ and 

z"{x) = y: (3) 

This summation is over the set Q^^ of rooted planar IPI graphs. 

This function, f{y), takes a rooted IPI blob and glues an arbitrary number of trees, 
each weighted with a factor of y, on to the blob. Note that each time we glue a tree on 
to a link in the blob, we pick up a factor of x, for the new vertex that is created, a factor 
of two because we can hang the tree in one of two directions and a factor of (1 + e^^)^ to 
take account of the 2^ different ways the spins on the blob and those on the tree's root 
blob can differ — this accounts for the factor of A multiplying each y. When we sum over 
all possible rooted IPI blobs and all possible ways of gluing trees to links this gives the 
first term in (^). More specifically, the blob without any trees attached has a weight of 
Z^\x). Consider a rooted graph with vertices, it has L = ^(3A^ — 1) internal links 
and we can add an arbitrary number of trees to each link, giving a total contribution of 
+ Ay + (Ay)^ + ■ ■ for the graph. When we sum over all IPI graphs this gives a 
contribution of 

- Ay)-^ = Z" (x{l - Ay)-i) . (4) 

The second term in (||) comes from the special case in which the root blob only consists 
of a single vertex, in which case we get a factor of x for the vertex, (1 + e"'^^)^'^ for the 
possible ways the spins may differ and y^ for the two trees. 
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Figure 2: T = f{T) 




The partition function, T, satisfies T = f{T) and this is represented diagrammatically 
in fig. 1^; the shaded circle represents the sum over IPI graphs, Z^\x). For fixed p and /?, 
this equation gives T as a function of x. By finding the closest singularity to the origin, 
which occurs at Xc{/3,p), say, we can determine the free energy for the model. For what 
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follows it will be convenient to define x' = x{l — AT)~2. Now, there are essentially two 
different types of singularity. One occurs when x' equals x^J , the critical value of x' for 
Z^\x') and this will correspond to the magnetized phase. The other occurs when the 
graphs become tree-like and we shall come back to this case later. 



2.2 Xc in the magnetized phase 



Consider the first case, for which x' = xjJ and hence Z^^{x') = , then writing T{xc 
as Tc, 

(5) 



Tc = y/l - XcTc Zi' + -^T2 



4rc. 



and xf = Xc(l - XcTc) 2 , where Ac = 2xc{l + e Defining, 



Xr\3 



we get 



2xf /i^ + -{l-h^ 



However, Z^\x) is a known function, from Brezin et al [10| one can show that 
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Z"{x) = -t(1 - 3r), 



(6) 
(7) 

(8) 



where = r(l — 2r)^. This gives Tc = g, x^^ = \J~^ and Zl^ = Hence we can 

solve (|^) to get Xc = x^^h? with 



^ -id 



d-l + J{l-ld){l + d) 



(9) 



where we have defined d = 

+ This gives Xc in the magnetized phase — we shall 

prove that it is magnetized later. 
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2.3 Number of blobs and '-)str 



Next we will calculate (n)^, the average number of IPI blobs at distance r. Define 
Gr = f^^\vT), where the notation means f{f{f{...f{vT)))); this weights each blob at 
distance r with an extra factor of v. Note that vT can be regarded as the standard 
partition function, but with the root blob weighted by an extra factor of v. Applying the 
function / to this gives f{vT), where now each blob at distance one is weighted with the 
extra factor of v. Each application of / just moves the weights of v down the tree, by one 
unit of distance. Thus, 



{n)r 



dGr 



Gr dv 



1 dGr 



v=l 



T dv 



where we have used that Gr{v = l) = T. However, G 



1 

T 



Of 



dy 



dGr-i 



dl 
dy 



y=T 



y=Gr-i 



dv 



v=l 
- v=l 



and thus 





'Of 




1 ~ 


dy 


y=T_ 



since {n)^ = 1 (there is only one root blob). It will be convenient to define, 

df 



B{x,P,p) 



dy 



y=T 



Using the formula (0) for f{y), one can calculate B, 



B 



dy 



2 ./I - \T V ^ ' ) ^ \ 



where 



T 2 VI -AT 

M"{x) 



+ e 



-2/3 



2p, 



dZ 



II 



Z^\x) dx 



1 - 3t 



(10) 

(11) 
(12) 

(13) 

(14) 
(15) 



Note that M^^(xf ) 
this gives us. 

Thus at Xc, {n)j 



|. Evaluating at Xc will give Bc{l3,p). In the magnetized phase 



Br 



l + ^(rf-l) 



1 



(d-l) 



{BcY and the average total number of blobs is 
in) ^ 



1 - Br 



(16) 



(17) 



For large enough /?, Be is less than one and positive, that is, the number of blobs decreases 
exponentially with distance. As /? is reduced, Be increases (since A increases and this 
encourages the tree to branch), until Be = 1 at some critical value, (3*; this corresponds 
to the boundary of the tree-like region. Using T = f{T), one can show that, 



dT 
dx 



1 



1- B 



d^T^ + ^''^^'^ ( M"{x') + xdT ( M"ix') - 1 
xyl — AT ^ ^ 



(18) 
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For j3 > (3* , < Be < 1 and ^ is finite, whilst for /3 < ^ diverges as x ^ x, 



dx 



That is, Be = 1 throughout the entire tree-hke region and 'jstr > in this region, since 
^ ~ {xc — x)'"'"'^^ . One can easily calculate '^str in the various regions of the phase 
diagram. In the tree-like region and on the boundary, |g ~ hence 1 — i3 ~ (Ax)'''^*'', 
where Ax = Xc — x. However, we also have an expression dll) for B and this gives 1 — B ^ 
(Ax)^"''^^*'' inside the tree phase. Hence, 7str = \ here as expected. On the boundary 
I - B (Ax)^^^"'^'''"^ giving 7*^^ = \. In the magnetized phase ~ ^^^al''^'^ ~ 
(Ax')^2 (Ax)~2j giving ^ystr = —\ as expected. 

Equation ( p!6| ) applies in the magnetized phase and also on the boundary between the 
two phases. At /?*, i?c = 1 and this gives using (p^), 

/3* = -iln(2F-l). (19) 

Note that this gives an estimate for the location, in the full model, of the transition 
between the tree and magnetized phases; compare this with the estimated location of the 



transition from the tree-like to the unmagnetized (pure gravity) phase, given in |12]. 



2.4 Xc in the tree phase 

Next we will calculate Xc in the tree phase. At Xc we have from Tc = f{Tc), 

Tc = - AcTc Z"{x') + Xc (1 + e-^^fr^ (20) 
and from Be = I, using (|l4|), 

where x' = Xc(l— AcTc)~2 . Remembering that d = (^1 + e^^^^^ and defining X = l—X^Tc, 
(H) gives 

Z" (2,M" - 1) x'Xd + 2d^Tex'X2 = 1 

and (pO) gives 



22 



Thus, 



Moreover, 



and thus 



d^T'^x'x'i =Tc-X^Z". (23) 



Z" (sm" - 1) x'Xd + 1 - Y^^^" = 0- (24) 



1 -X 1 -X 

Tc = = 3- 25 

Ac 2x'X2d 



Z" (sm" - 1 ) X{1 - X)xd +1-X - Ax'Z"x'^d = 0. (26) 
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Also gives 

Ax'Z"X^d ={1-X)[2- d{l - X)] . 
Combining these last two equations gives 



Z"{3M^^ - l)x' - 1 



Xd + d-l = 0. 



Using (|) and (|l5|), 



and substituting into p7|) gives 



and thus Xc = x X2 is given by 



X 



1 r 



d-1 
d{l - 2t) 



l + 2((i- 1) 



1 . 3 



Xc = -d 2 — 1 



(27) 

(28) 
(29) 

(30) 
(31) 



in the tree phase. From equations (31) and (^) one can easily show that there is a third 
order phase transition with a finite discontinuity, that is, the critical exponent a = —1. 



2.5 Number of vertices 



One can also calculate the average number of vertices at distance r. Define = 
j{r) (^f(^T)\xz), where the notation f{y)\xz means the function f{y), but with x replaced 
everywhere by xz. Thus f{T)\xz is just a tree, but with each vertex in the root blob 
weighted by an extra factor of z. Applying /, r times, just pushes these weights down 
the tree to a distance r, so that in all the vertices at distance r have an extra weight 
of z. Thus the average number of vertices at distance r is given by 



' z 


1 




[G'r dz \ 


z=i~T 


_ dy 



dz 



Z = l 



(32) 
(33) 



and this is essentially just the (geometric) two-point function (see section ^). The average 
total number of vertices is 



where {N)q is the average number of vertices in the root blob. This is given by 

(N)^ = B + M"{x') - xd^T (m"{x') + 1) . 
In the magnetized phase this gives, at Xc, 

Wo 



1 



{d - If 



(34) 



(35) 



(36) 
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and in the tree phase, 

Note that d = (1 + e~^^)^, with 1 < d < 2 in the magnetized phase and c? > 2 in the 
tree phase. The formulae also show that the average number of vertices in a blob is very 
small, for any values of f3 and p; in fact {N)q = 2 at /?*. 



3 Critical exponents from the scaling relations 

Prom equation (^) we see that there is an exponential decay of the number of blobs with 
distance. One can define a mass, m, in the model through, 

(n)^ = B'- = e-"^^ (38) 

The mass is a function of /3, p and x. At Xc, as the critical line is approached from 
the magnetized phase, the mass vanishes (i.e. there is a correlation length equal to 1/m, 
which diverges). Note that at Xc, m = throughout the whole tree phase. Let us use 
this definition of the mass and various scaling relations to calculate the critical exponents. 
Later we will calculate the magnetic exponents directly from the model, for the case p = 1, 
gaining the same results for (3m and 6m', this will show that our definition of the mass is 
reasonable and that the scaling relations hold for this toy model. 

First let us consider the geometric exponents; note that we have already calculated ^str 
in section |2.3| . The exponent u is defined by m ~ {Ax^ for Ax (where Ax = Xc — x), 
but m = —InB so that m ~ 1 — i3 ~ (Ax)^"'^'' and hence u = ^str] in the tree phase 
= Istr = \ and on the boundary v* = 7*^^ = |. Since there are no power law corrections 
to ( |38|) we have r/ = 1; note that at small Ax, for 1 ^ r ^ 1/m, one might expect []TT| 
that (n)^ ~ r^~'^ . These sets of exponents satisfy Fisher's scaling relation ^str = ^^(2 — t]). 



In reference |11] it is shown that if the two-point function has associated with it a 
vanishing mass, then for a suitable definition of the Hausdorff dimension, udn = 1- This 
gives that dn = 2 in the tree phase, as we might expect for branched polymers and d*jj = 3 
on the critical line. 

Consider now the magnetic exponents, which we will write with a subscript m to 
avoid confusion. Evaluating at Xc and letting /3 — > /?* from the magnetized phase, we 
have, defining A/3 = (3 — {3* , m ^ Af3 from ( |l6|) and hence Vm = ^- If we take rjm = I 
(we shall see later that the spin-spin correlation function has no power law corrections 
either), then using 7^ = fm(2 — rjm) gives 7m = 1. Applying the other scaling relations 
2 - Q = i^md*fj, (3mSm = Pm + 7m and a + 2Pm + 7m = 2, yields a = -1, /3m = 1 and 
6m = 2. The various exponents are listed in table |l]. 
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Table 1: Critical exponents 



Phase 


Istr 


U 


rj 




Magnetized 


1 
2 




(1) 




Critical 


1 

3 


1 

3 


1 


3 


Tree 


1 
2 


1 
2 


1 


2 



a 




7m 


<5m 






-1 


1 


1 


2 


1 


1 



Before proceeding to check the magnetic exponents explicitly, we shall define dn and 
calculate it directly from the model. 



4 Hausdorff dimension 



As before in |11| we define the Hausdorff dimension, dn, in terms of the two-point function 

T2{r)= ^"^G, (39) 

Geg|(r) 

where Q2{f) is the set of planar (jy^ graphs, with two marked points that are separated 
by a distance r. One of the marked points will be taken to be the vertex, in the root 
blob, which is connected to the external leg. Wq is just the usual weight for the spin 
configurations, which appears in (|^). Then du is defined, in the continuum limit, by 



oo. 



where 



N{r) = 



T2{r) 



N Wg 



m{Ax)r = const., 

a(lnT2(r)) 
dx 



(40) 



(41) 



That is, by tuning x to Xc we are taking a continuum limit; however this definition of dn 
only really makes sense in the tree phase and on the boundary, where the mass vanishes 

Now T2(r) is just 



T2{r) = Y x^WcNr = T{N)^ = TB^{N)^, 



(42) 



where is the number of vertices at distance r for graph G; 

N{r) = X 

As r — > oo, 



dM" dT 
N{r) ~ r ^ . ~ r 



IdT r dB 1 d{N)Q 
f~dx ^ B~dx ^ {N)q dx 

(2;f -x')~^(Ax)-^=' 



dx' dx 



(43) 



(44) 
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Since we are taking the continuum limit with mr fixed, r ~ m ^ ~ (Ax) , 

iV(r)~r2(xf (45) 

In the tree phase N{r) ~ giving du = 2 and on the critical line N{r) ~ r^(Ax')~2 ~ 
r^(Ax)~3 ~ r^, so that = 3, as expected. Unfortunately, due to the way distances and 
the Hausdorff dimension have been defined in this model, dn is not well-defined in the 
magnetized phase. In reference [ p]] it is shown, using the geodesic distance, that dn = 4 
for pure gravity (also the geometric exponents are v = 1/4 and r] = 4) and this may well 
be true in the whole of the magnetized phase. It is certainly the correct value at p = 0, 
where the toy model reduces to a pure gravity model. 

Thus our definition of distance appears to correctly capture the behaviour of the model 
in the tree phase and on the critical line, where 'jstr > 0, but not within the magnetized 
phase. 



5 Spin-spin correlation function 

In this section we will calculate (SS)^. the average total spin at distance r; the summation 
is over all vertices at distance r and all spin flavours on those vertices. Since the spins 
in the root blob are fixed to be +1, this is essentially a spin-spin correlation function, for 
spins separated by a distance r. To calculate this we will add a magnetic field for vertices 
at distance r. First we will solve the model with p = 1 and then the general case. Let us 

define 

U{y) = - Xye" (^xe^ (l - Aye^)"^) + x(l + e-'f')Ve'' , (46) 

with A = 2x(l -|- e~^^). This is just f{y), but with x replaced by xe^ . The function takes 
a number of trees with weights y and glues them on to a blob, where each vertex of the 
blob has a spin of -|-1 on it and there is a magnetic field +H applied on these spins. The 
function f~{y) is defined in the same way, but with H replaced by —H throughout. 

Let be a tree with positive spins in the root blob and a magnetic field of +H 
applied at distance r from the root. Similarly T~ has spins of —1 in the root blob. Then 
T+ = U{T) and = f.{T). Define, 

F{y) = VI - 2xy Z" {x (1 - 2xy)-i) + xy^. (47) 

Then 

r+ = F(T;ti + e-2^T-_i) (48) 

T- = F(T-_i + e-2/^r+i), (49) 
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that is, to make a graph with spins of +1 on the root blob and a magnetic field at distance 
r, one takes a number of trees with magnetic fields at distance r — 1 and glues them on 
to a blob with positive spins, picking up a factor of e~^^ if the spins on the root of the 
tree are negative. Now, 



1 dT+ 
dH 



1 dT+ 



H=0 



T dH 



(50) 



H=0 



the superscript on {T,S)^ denotes the fact that the vertices in the root blob have a single 
spin each, which is fixed to be +1. 



dH 



dH 



+ e 



dH 



dF 
dy 



(E5),+ = {ES):_, + e-'^{ES);_, 



y=T+_,+e-^0T-_, 
dF 

dy 2/=(i+e-2/3)T 



(51) 
(52) 

(53) 



but {^S);_^ = and F{y) = f{2xy/X) so that 

dF _ 2xB 

dy 2^=(i+e-2/3)r A 

Thus defining t = tanh /?, 

= Bt = iBtni:S)^ = {BtYiN},. (54) 

Let us consider the general case of p spin flavours for which X = 2x 

T++-+ = F(r+_:i+ + (r+_:i+- + t+_:;+-+ + •••) + 

e-^^ (t+_;i+-- + ...)+... + e-^'^^'T-Jr) (55) 

2xB 



(S5) 



++•••+ 

r 



A 



[(S5),f 1+ + e-'^ ((SS),f i+- + •••) + ••• + e-'^P{^S);:r] . (56) 



However one can easily show that. 



1 



(ES);- = -(p-2n)(S-S) 



++•■ 

r 



(57) 



where there are p—n plus signs and n minus signs in the superscript on the left-hand side. 
Thus, 



++•••+ 

r 



= (i?i)(S5)++-+ = {Btyp{N),. 



(58) 
(59) 
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Now defining the magnetization at distance r by 

Mr . (.0) 

we have that A^,. = V . It should be noted that if one uses the exponential decay of 
{^S)^ ^ or M-r to define a spin-spin correlation length, then this quantity will not 
diverge at the phase transition; this behaviour should be contrasted with that of the 
geometric correlation length, which does diverge. 

If we define the total magnetization M , for the case in which all the spins in the root 
blob are positive, by 

where (SS*)^ ^ is the average of the total spin, which is a sum over all vertices and all 
spin flavours, then we have, at Xc, 

^ = («^' 

Note that M = 1 for /? = oo, and that = throughout the tree phase, showing that 
this phase is unmagnetized. The other phase has < Be < 1 and hence < < 1; 
this is the magnetized phase as claimed earlier. Also at Xc, near 1 — Be ^ A/?, where 
A/3 = /? — /?*, so that M ~ A/3 and thus Pm = as calculated from the scaling relations. 
In fact for finite p, 

M = fAp + 0{{Apf), (63) 

so that the coefficient in front of A/3 is non-zero in general (note that p = would give 
(3* = -oo). 



6 Toy model, with a magnetic field 

In the previous section we defined the magnetization in the grand canonical ensemble by 
(|6l|). A different definition is possible namely, 

A^,e= lim -^{^S)j„ (64) 

N-*oo pJM 

where now we are working in the canonical ensemble, that is, we are using the set, QliN) 
of rooted A^-vertex graphs. 

In this section we will add a magnetic field, H, to the toy model. However, we shall 
only consider the case p = 1, as p > 1 appears to be more difficult to solve. It will be 
shown that the magnetic exponents /3m and 5m are the same for both definitions of the 
magnetization and agree with the results from the scaling relations. 
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Rather than use the transfer function f{y), which glues trees on to a IPI blob, we 
shall use a different transfer function /(y), that glues trees on to a domain; this is more 
convenient for the p = 1 model with a magnetic field. For H = 0, it will give a sum over 
the same set of graphs that we had previously, with the same weights, but we will no 
longer be able to keep track of distances within the graphs. 

Consider first the case H = 0, the function /(y) is given by 

m = ^Jf^^ (1 - 2xKri) + xK\ (65) 



where 



with = 4x6 and 



K 

2x 



^ (l - v/r^) , (66) 



Z\x) = Y: (67) 



The function defined in (65) takes a rooted (p^ graph, which will form the core of a 
domain, and glues on trees weighted with factors of K. As before one picks up a factor 
of 2x for each tree glued on. In this case however we allow trees to be glued on to the 
root link, which changes the power of the factor in front of compared with that in (Q). 
This is necessary to make sure that we correctly sum over all possible domain structures. 

The factor K is the solution of 

K = e-^^y + xK^. (68) 

This generates tree graphs whose vertices are weighted with x and which have a factor of 
e~^^y at the end of each branch. The domain corresponds to the core (f>^ graph plus all 
the vertices in K. The ends of the branches in K, which we have weighted with e~^^y, 
correspond to the domain boundaries; setting y = T will give us a set of domains glued 
together in a tree- like fashion. 

Thus we have 

Ry) = ^—-T ix (1 - /.y)-t) + ^ (l - v/r^)' (69) 

and the grand canonical partition function ([l|) is the solution of T = f(T). This is shown 
diagrammatically in fig. |3|. 

Note that this is just a different way of formulating the same model that we solved 
previously. One can easily use T = f{T) to calculate Xc in the two phases, gaining the 
same results as before (to solve the p-flavour case one just uses fi = 4x((i — 1)). The 
advantage of this new formulation is that as we use / to move down the tree from one 
domain to the next, the sign of the spins alternates (at least for p = 1). This simplifies 
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Figure 3: T = f{T) 




+ 



+ 



K + 




-2 P 



+ 




greatly the task of adding a magnetic field. For the original version of the model it is 
quite difficult to keep track of the spins on the different blobs. 

Now let us add a magnetic field of +H, the transfer function which makes a domain 
with spins of +1 is 

2 



(1 



1 



f+iy) = - ^^Z^(xe^(l-Mye^)"^)+^(l-Vl-^2/e^) • (70) 



Axe 



Each spin in the domain is given an extra weight of . It will be convenient to define 



xe 



H 



{l — fiye^^ ^ . Similarly f-{y), which makes a domain with negative spins 
is (70), but with H replaced everywhere by —H. Let be the partition function for 
graphs with positive spins in the root domain and T" that for those with a negative root 
domain. Then T+ = /+(r~) and T~ = f-{T^). These equations determine T^{x,/3,H). 
For if > in the magnetized phase, the critical value of x is determined by x'_^ = x^, 
where x^ is the critical value for Z^{x'_^_). Defining 



XrG 



(71) 



we have 



and 



(72) 



(73) 



Note that x[ = 1/(2.33), = 33 (l - ^) and = l/(\/3- |) - 1 from @. However 
we still have T~ = f^{T^), and using this gives an equation determining h. Putting 



X = 1 



IJ.ce 



-H 



(74) 
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and x' = Xre X 4 then we have 



1 + X + 2VX [2x'Z^{x') - I 



This gives h{(3,H) and hence Xc{(3,H). Now, 

1 dxr 



At ^ we get by differentiating 

MceiP,H^O+) = 1 



Xc dH 



3 dh 
^ ~ hdH' 



3 (2 - - ^(h) 



-2 



(3 - h'^VS - (2e2/3 + e-2/3 + 3) 



with 



2(1+6-2/3) 

note that at = 0, /i2 = 2/i2/a/3, see equation (|). At /3 close to /3* = 0, 

7W,e(A/3, H^0+)= 2A/3 + 0((A/?)2). 



(75) 



(76) 



(77) 



(78) 



(79) 



Hence, Pm = 1 as previously, but note that Aice is not equal to M.. Now consider dm, 

1 

defined by A4ce(/?*, A//) ~ (A//)''™. From (ffSD, one can show that Mce is a certain 
function of h, e''^'^, e^^ , x' , Z^{x'), M^(x') and X, where 

x' dZ\x') 



M\x') 



ZHx' 



dx' 



(80) 



At (3* and small /S.H, barring accidental cancellations, then all these functions are a 
constant plus terms of order /S.H (or smaller), except for M^(x'), 

M\x') = Mi+(xi-x')^ + ■■■ . (81) 

Since x' = + 0(Ai?) then M\x') = + 0{{AH)^), giving 6m = 2. Note that the 
value of 6m depends crucially on the fact that ^str = — ^ in the magnetized phase. 

Suppose that we next try to calculate 7^, defined by 



(82) 



then we find that x ~ ^-iid that this contains terms such as ^ g^rf ^ ; which diverge 
as we take H to zero. Thus jm seems not to be well-defined, even though the scaling 
relations give = 1- 

Consider now the magnetization in the grand canonical ensemble, M, defined by (pl|), 
f3m has already been calculated and this just leaves 6m- We see that Ai ^ at P*, due 
to the divergence of (N), 



MiP*,AH) 



(N) 



(83) 



/3* 
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For graphs with a positive root domain, at Xc, 

T+ dx 

where T+ = /+(/_(T+)). This gives 



dx 



V-f 




df- 




dy 


T- 


dy 


T+- 



Thus 



MiP*,AH) ~ 1- 



dx 



+ 



Of, 



dy 



df- 



dx 



T+ 



dy 



df- 



dy 



(84) 



(85) 



(86) 



and again this is a function of various variables all of which are a constant plus 0{AH), 
except that M^(x') = + 0((Aif)i), giving 5m = 2. 

This completes our analysis of the toy model; we have derived all the critical exponents 
from the model (except for 7^) and shown that all the usual scaling relations hold. Now 
we shall turn our attention to the full model of p-Ising spins coupled to 2d gravity, and 
will find that many of the considerations in the previous sections apply for the transition 
between the tree and magnetized phases in this model. 



7 Full model 



7.1 Definition 



Consider the full model where we have p independent Ising spins on each vertex and are 
summing over all rooted planar (p^ graphs and all spin configurations. Then the partition 
function is, with t = tanh/?, 

p 



Z\x,(3) 



N 



1 

2^ 



(87) 



{5} <ij> 

where we have divided out various factors of two and cosh/? in order to simplify the 
formulae later on. N is the number of vertices in graph G. Note that there is no vertex 
and are no spins at the end of the root, so that the product over links does not include 
the root link. In fact we shall consider a generalization of this model that has p coupling 
constants. Expanding out the above equation gives, with the extra coupling constants. 



1 



{S} <ij> 



p-1 



Z^(x,{A})= ^ x^^^ n l + XiEsrS]' + X2Yl E SrS-S^S^ + ... (88) 



a=l /3=a+l 



+Ap n ^T^: 
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where Sf is the flavour a spin on vertex i, the second summation is over aU spin config- 
urations and {A} is the set of couphng constants {Ai, A2, • " " ) Ap}, with < Aj < 1 for 
ah j. Of course by setting Aj = P one recovers the ordinary model (^)- Alternatively by 
setting Aj = for j > 1 we recover the 0(n) models studied in |l^ with n = p. To 
save writing we define T = Z\x, {A}). We will show that T satisfies T = f(T) where 

f{y) = y^T^ Z"{x\ {A'}) + xy^ (89) 

with 

3 



x' = x(l-2xy)-2, A; =A,- ; (90) 



(1 - 2xy) 
(1 - 2xy\,) ' 

Z^\x\ {A'}) is defined as in (|88|), but using the set of rooted IPI graphs, Ql^ . In a similar 
fashion to (|^) the function f[y) takes a IPI graph and glues trees on to it, but this time 
as well as x being changed to x' , the coupling constants are also renormalized. Note that 
a similar renormalization is studied in |p. 

To justify the above equations let us first consider the case p = 1. For each link we 
have a factor such as (1 + XiSiSj) and when we multiply these factors together, the sum 
over spins will cause any terms containing odd powers of a given spin to vanish. The 
only non-vanishing terms correspond to sets of closed loops; we shall refer to such a non- 
vanishing term as a loop configuration. Thus for each graph G, we sum over all ways of 
drawing sets of closed non-intersecting, non-backtracking loops on the graph (call this set 
of loop configurations Cq)- So, 

Z'ix,Xi)= E E ^1' (91) 

where N is the number of vertices in graph G, and / is the number of links making up 
the loops in the loop configuration C. 

As before, we will make the graphs G (where G £ Q^) by gluing together IPI blobs 
(fig. H). Note that because of the tree-like structure, any given closed loop is wholly 
contained within a single blob. This means that the partition function for a given graph 
G, factorizes into a product of contributions from the individual blobs. To calculate the 
transfer function, /(y), we take a rooted IPI blob, summing over all graphs and loop 
configurations, as well as over all ways of attaching trees to the blob. If no trees were 
attached, one would just have Z^\x,\i), which is 

Z'\xAi) = E ^''wE n [1 + Ai5.5,] = E E Ai 1^"'' (92) 
G&g? {5} <ij> Gegi' c&Ca 

where L is the number of internal links in graph G (i.e. excluding the root link); that is, 
L = ^(3A^ — 1). Attaching a tree weighted with y, gives a factor of 2xy (there is an extra 
vertex and there are two ways of hanging the tree off the link). For a given G and £, 
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each link either contributes Ai if it is part of a closed loop or 1 if it is not. Attaching an 
arbitrary number of trees to an internal link that contributed 1 changes this contribution 
to {l — 2xy)~^; however for those that contributed Ai we get Ai(l — 2xyAi)~^, since adding 
a vertex on a link which was part of a closed loop, increases the length of that loop by 
one and hence gives an extra factor of Ai. Thus summing over all ways of attaching trees 
causes the change: 



1 



Xi{l-2xy) \ 



and hence 

Z"{x, Ai) — > ^1 - 2xy Z"{x', X[). (94) 

The extra term in (^) comes from the case in which the blob is just a single vertex. 

The general case for p > 1 follows with only minor modifications to the argument. 
Each closed loop is labelled with a spin flavour (an integer a, with 1 < a < p) and 
loops with different spin flavours can intersect. Thus for a given graph G and loop 
conflguration C, a given link will have j spin flavours running through it (0 < j < p), 
giving a corresponding factor of Xj. Adding an arbitrary number of trees to this link 
changes the contribution to Aj(l — 2xyXj)^^. Hence, we get (|8^) and (|90|) for the general 
case of a rooted IPI blob with trees, each weighted by y, hanging off it. Setting y = T, 
we recover the partition function for graphs with G G Q^, that is, T = fiT); f[y) is the 
transfer function for the full model. 



7.2 Exponential decay of blobs 

As before, we define distance to be the geodesic distance, but count only the links con- 
necting IPI blobs. Again one has Gr = f^^'HvT), giving the partition function for a 
model where all blobs at distance r are weighted by v. The average number of IPI blobs 
at distance r is 



dGr 



Gr dv 



1 dGr 



v=l 



T dv 



(95) 



v=l 



Since Gr = f{Gr-i), 



1 df 


dGr-l 




dl 


T dy 


y=T 9v 


v=l 


dy 



in) 



r-l' 



(96) 



so that (n)^, = with B = ^ ^ . Thus we still have an exponential decay of the average 
number of blobs with distance and 



B 



xT 



l-2xT 



1 - 3xT + 



(1 - xT) 



BZ" ^ dZ" 



dx' 



dX!^ 



(97) 



y=T 
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where Z^^ = Z^\x' , {A'}). One can also show that 

1 



( df{y) \ 
[l-B)\ dx y|A} 



x(l - B) 

{Xc — x)~^''^ 



y=T 

B-xT+{l- xT)M"{x', {A'}) 



where 



M"ix',{X'}) 



dZ 



ir 



Z^^ \ dx' 



(98) 

(99) 
(100) 

(101) 



{A'} 



Thus as in the previous case, at Xc, B^ = 1 in the tree phase and on its boundary, where 



Istr > 0, and we have < Be < I in the magnetized phase. Using equations (|97|) and (|99| ) 
one can show, barring accidental cancellations, that jstr = ^ the tree phase and that 
on the boundary, 7*^^ = 7^^/ (7^^ — 1), which is just Durhuus' formula where 7^^({A'}), 
which is negative, is the value of the string susceptibility for Z^\x\ {A'}), that is, 



d^z 



11 ' 



dx''^ 



(4({A'})-x')'(^+^"«"'»^ 



(102) 



{A'} 



7.3 The spin-spin correlation function 

Next we consider the spin-spin correlation function. Rather than fixing all the spins in 
the root blob, as we did for the toy model, it is more convenient to add a set of p external 
spins on the root, which are all fixed to be +1. The partition function is as in (p8|), but 
with an extra factor for the link connecting the external spins to the root blob, which we 
shall denote Wext (we will not add an extra factor of x). Thus, 



Z\x,{\}) 



Gag', {S} 



W, 



ext n 

<'J> 



(103) 



where the weight Wij for link < ij > is given by the square bracket in equation ( pSj ) . Note 
that when the sum over spins is performed, the factor Wext just gives a contribution of 
one, so that the function Z^ is unchanged by adding the external spins. The average total 
spin at distance r is 



2,1 ^ oNv 



2Np 



Wextl n WiAY^s^ 

.<ij> I k,a 



(104) 



Geg^ {S} 

where the last summation is over all spins, k, at distance r and all spin flavours, a. 

Consider Z^ (T^S)^ , the only non-zero contributions come from the configurations in 
which there is a flavour a line from the external spins to vertex k, in addition to the usual 
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sets of closed loops. Define Gr by 



Gagl {S} 



.<ij> 



pGr, 



(105) 



Then Gr is a sum over all rooted graphs and all loop configurations that have a flavour 1 
line from the external spins to a vertex k at distance r; the location of the vertex k is also 
summed over. However, we can write Gr in terms of Gr-i- This is shown diagrammatically 
in fig. 1^, where links included in the flavour 1 line are drawn thicker than those which are 
not included. The symbol on the left-hand side is used to represent G^, it has a flavour 1 
line of length r, note that we are still measuring distances from the root blob (not from 
the external spins). 

Figure 4: Formula for Gr 



r-l 




Gr 



(106) 



That is, 

Gr = hXixT + Zi\x, {A'})Ai(l - 2xT) 

where ^2^(x, {A}), which is drawn as a shaded circle with a thick line passing through it, 
is the partition function for IPX graphs with two distinct legs and a flavour 1 line running 
through (it does not include factors of Ai on the two legs); x' and {A'} are given by (|90| ) 
evaluated at y = T. Define D to be equal to the square bracket above, then Gr = D^Gq. 
So that we have an exponential decay with distance. 



The total average spin is 



1 



l-D 



The magnetization in the grand canonical ensemble is given by 

_ /I - s\ .. n-B 



M 



p{N) p{N), 



1 - D 



Mr 



l-D 



(107) 



(108) 



(109) 



where A^o is the magnetization of the root blob (0 < Mo < 1). We have used the fact 
that the average number of vertices at distance r is (A^)^ = B^{N)q, so that (A^) = 
(A^)q/(1 — B). This is easily shown by defining G'^. = f^'^\f{T)\xz) and following the 
derivation of (H). Note that, at Xc, in the tree phase Be = 1, so that M = throughout 
this phase as one might expect. 
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In the next section it is assumed that D < 1, that is, that the spin-spin correlation 
length does not diverge; the easiest way to justify this would be to show that D < B, 
as we already know that B < 1. The following calculation will make this assumption 



plausible. From (107) one clearly has D < B, {at least provided that (SS')q 0), since 
(SS"),^ < p{N)j.. To improve this inequality we need an equation for B which is simpler 
than (|97|). Define G" to be the partition function for rooted graphs, with a marked vertex 
at distance r — it is a sum over all graphs, G ^ Ql, all loop configurations and all ways 
of marking a vertex at distance r; that is, G'l = T{N).^. Then following the derivation of 
(llO£| ) we have 

G;' = hxT + Zl\x\ {A'})(1 - 2xT)l G;'_i = BG",^^, (110) 

where Z^^ is the partition function for IPX graphs, with two distinct legs — it is a sum 
over all such graphs and all loop configurations. Note, Z^^ differs from Z2 in that there 
is no flavour 1 line running through it. Thus, 

B = 2xT + Zl\x\{>!]){l-2xT) (111) 



D = \i 



2xT + Zi\x\{\'}){l-2xT) . (112) 



Now, if one could show that 2^/^(x', {A'}) > Z2 {x' ,{\'}), then this would imply that 
D < XiB. 

Consider Z^_{x' , {X'}), the partition function for IPX graphs with two distinct legs, 
in which the flavour 1 spins on the two vertices connected to the legs are held fixed 
at +1 and —1 respectively; all other spins are summed over. Then one can easily see 
that Z^_{x' , {X'}) = Zl^{x',{X'}) — Z2 {x' ,{X'}), since the only non-vanishing terms 
correspond either to sets of closed loops, or configurations in which there is also a flavour 1 
line running between the two fixed spins; this last set of terms is multiplied by the product 
of the fixed spins, namely —1. For small enough values of {A'}, Z^_ [x' ,{X'}) is manifestly 
positive; the weights Wij for each link will be positive for any spin configuration (see 
equation (|8^)), and hence Z!^_ is just the sum of positive terms. At X'j = 1 for all 
j, the system is completely magnetized and hence Z^_ = 0. One would expect that 
Z^_{x' , {X'}) > for other values of {A'}. Xn which case D < XiB, with D < B in 
general and D = B only when Aj = A^ = 1 for all j. Thus in general D < 1 and the spin- 
spin correlation length does not diverge. Note that our definition of distance is only good 
in the tree-like phase and on its boundary, and hence this result does not say anything 
about the possible divergence of the proper spin-spin correlation length (defined using 
the geodesic distance) at the magnetization transition, between the U and M phases (see 
figure 1). 

Figure ^ shows a tentative phase diagram for the standard model defined by (^) ; T 
is the (unmagnetized) tree-like or branched polymer phase, U the unmagnetized phase 
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and M the magnetized phase. See 13] for a discussion of various possible alternative 



phase diagrams. Each line of this phase diagram for constant p (where p is a non- negative 
integer) corresponds to a line through the p-dimensional phase space of the generalized 
model. Within the U and M phases the model behaves in a similar fashion to the pure 
gravity case and it seems probable that jstr = throughout both these phases. The 
transition between the magnetized (M) and unmagnetized (U) phases is caused by A4q 
vanishing as the coupling constants {A} are varied. From the KPZ results, we expect to 
have — ^ < 7str — along the corresponding critical line, at least for p < 2. The tree-like 
phase has ^str = | throughout, the generic value for branched polymers. 




7.4 Critical exponents 

Now consider the transition from the magnetized to the tree phase, which is caused by 
Be ^ 1. The behaviour of the full model at this transition is very similar to that of the 
toy model. Following the analysis in section |^, we define a mass hy m = —InB. As 
before there is an exponential decay, so we will take rj = 1, which gives jstr = v provided 
that the scaling relation ^atr = ^{2 — ij) holds. Consider the magnetic exponents, at 
Xc and fixed integer p we vary the {A} in order to approach the phase transition; the 
variable AA will be used to parameterize this, with AA = at the transition. Then by 
definition m ~ (AA)'^™, but m = — In Be so that Be ^ I — (AA)'''". Now /J^ is defined by 
M ~ (AA)^™, but 7W ~ 1 - 5c from (|09|), thus pm — Ufa] note, we have assumed that we 
are not close to any point at which A^q vanishes {P in fig. |5| is such a point). It has also 
been assumed that Dc 7^ 1 at the transition, i.e. that the spin-spin correlation length does 
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not diverge; this was justified in the previous section. If we assume that ^rn — ^m(2 — ^m) 
with r/m = 1, then we have 7m = i^m = Pm- Assuming the scaling relations 2 — a = I'md*^ 
and a + 2 (3m + 7m = 2, gives d*jj = 3. The calculation in [11|, which gave that udn = 1, 
still applies, so that u* = ^ and hence 7*j^ = |. This is encouraging as the analysis in ||8| 
leads one to suppose that 7sir — 3 

provided ^str = ~| i^i the magnetized phase. Also, 
applying the relation (3mSm = Pm + 7m gives (5^ = 2. 

Now B is given by ( |97[ ) and in the magnetized phase at f ({A'}). Provided 

that we are away from points at which A^o vanishes, Z^^{x',{X'}) should be analytic 
and Taylor expandable about {A'*} (the value of {A'} at the M to T phase transition). 
Assuming that Xc({A})Tc({A}) is expandable in powers of AA as far as AA, then barring 
accidental cancellations Be ^ 1 — A A and hence Vm = Pm = 7m = 1 and a = — 1. 

Thus on the assumption of various scaling relations and that our definitions of m 
and M are adequate, and barring problems such as accidental cancellations, the set of 
exponents given in table |l] should also apply for the full model for the M to T transition, 
away from point P. Again it is not clear what the value of dn is, within the M phase; 
however one might guess that dn = 4, the pure gravity value [0, 16, throughout the 
U and M phases (in which case it is probable that v = 1/A and r] = 4 here). Note that 
the mechanism by which the U to T transition occurs is essentially the same and so we 
would expect to have the same set of exponents, although the magnetic exponents are not 
relevant since = in both the U and T phases. Perhaps one should note that a low- 
temperature expansion matrix model calculation in the limit c — > cxd gives that 7*^^ = | 
at the transition between the tree and magnetized phases, and that the truncated model 



studied in ||12[, which is supposed to approximate the U to T transition, has 7*^^ = 3 
and a = — 1, agreeing with our predictions. In addition the q-state Potts model at large 
q j5[ , which is equivalent to the multiple Ising model in the limit q = p = 00 and may 
be related to it for finite q, has a branched polymer region (with jstr = ^) separating two 
pure gravity regions (with ^str — ^^'^ has 7*^,. = ^ and a = —1 at both transitions. 

At point P, the above argument fails since Aio — > and hence it is possible to have 
Pm 7^ i^m; some of the other assumptions may break down here too. Thus at this point 
one might expect to get a different set of exponents. In any case the calculations in this 
section only apply for non- negative integer p and it is not entirely clear whether p* , the 
value of p at point P, is an integer. It is tempting to suppose that p* = 2, in which case 
one could understand the breakdown of the KPZ formula in terms of the appearance of 
the branched polymer phase, however Monte Carlo simulations do not seem to support 
this hypothesis 0]. 
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8 Conclusion 



In this article we have studied correlation functions in two-dimensional quantum gravity 
coupled to Ising spins. For the toy model two approximations are used: (i) the spin 
configurations allowed are the ones dominant in the low temperature expansion of multiple 
Ising spins on dynamically triangulated surfaces and (ii) the distance between two spins 
is identified with the number of links separating IPI subgraphs along a path connecting 
the spins. In the full model only the second approximation is made. We do not expect 
this last approximation to be important for values of /3 where jstriP) > 0, i.e. for values 
of P where the average number of vertices of a generic t/)^ graph diverges for xyxciP), 
since the average number of vertices in an IPI part of a generic 0^ graph in our ensemble 
is very small. However, if jstriP) < the average number of vertices in a generic (p^ 
graph is itself small, even at the critical point Xc(/3) and our definition of distance can not 
be used to extract in a reliable way the fractal properties of the ensemble of (p^ graphs. 
Approximation (i) has been shown to be exact in the c — > oo limit [Q] and numerical 
simulations |^ suggest that it is an excellent approximation except for the smallest values 
of P if many Ising spins are coupled to two-dimensional gravity. For moderate values of c 
it is difficult to judge if approximation (i) is reliable all the way down to f3* . For c = 1/2 
and c = 1 (i) not a good approximation for (3 £ [0,/3*]. 

The toy model allows us to analyze the two-point function (the puncture-puncture 
correlator) and the spin-spin correlator as a function of the distance r. For /3 < /?* we 
have jstr > and the fractal structure is reliably extracted from the two-point function. 
We find that dn = 2 for jS < (3* , while it jumps to dj^ = 3 at /?*. Whilst we can not 
calculate dn for f3 > (3*, one might expect it to equal four, the pure gravity value |11]. The 
fact that dn is different from 4 for /? < /?* is an indication of the very strong interaction 
between gravity and matter for c > 1. 

A further indication of the strong link between geometry and matter configurations 
present in the model is found in the scaling relations for the magnetic exponents. They 
are found from the behaviour of the two-point function {n)r in the region (3 > [3* . The 
exponential decay of (n)^ at Xc{(3) as a function of the distance r defines a length scale 

m ~ j^. (113) 

which diverges for (3 ^ (3* and allows us to define the magnetic exponents using the 
standard hyper-scaling relations (note that ^{(3) = oo for /3 < /?*). The same exponents 
are found, for the toy model, by a direct calculation in the canonical and grand canonical 
ensembles. Strictly speaking, it seems more natural to define magnetic scaling properties 
from one of the spin-spin correlators, i.e. either (SS)^"' or Air, however the corre- 
sponding correlation length does not diverge at This is an indication of the geometric 
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nature of the transition. 

For the fuh model, a summation over all spin configurations is performed and the only 
approximation made is in the measure of distance that is used. As before there is an 
exponential decay of both the two-point function and the spin-spin correlation function, 
and one can define, in the magnetized phase, a geometric correlation length, which diverges 



as the tree phase is approached. Based on the analysis in section [7^ we do not expect the 
spin-spin correlation length, at least as we have defined it, to diverge at the magnetized 
to tree-like transition. Using the geometric length scale and the scaling relations, which 
have been shown to hold for the toy model, all the critical exponents are calculated for the 
branched polymer phase and its boundary, on the basis of a small number of assumptions. 
The exponents are the same as those for the toy model (see table |l]) and again dfj = 2 
within the tree phase and dn is equal to three on the critical line, showing the strong 
interaction between the matter fields and the geometry. 

Some interesting questions remain to be answered, such as what happens outside the 
branched polymer phase, in particular what the values of the geometric exponents and 
dn are, and whether there is a diverging length associated with the spin-spin correlation 
function at the magnetization transition (between the U and M phases). Unfortunately 
our measure of distance is inadequate in this region and it will require further work using 
probably the full geodesic distance in order to determine this. It would also be useful 



to locate point P on the phase diagram, as it has been suggested Jl^, y] that p* > 2, in 
which case there is an intermediate region, for 2 < p < p* , between the KPZ regime and 
the branched polymer phase. 

Finally, it should perhaps be noted that most of the considerations in this paper also 
apply to the 0(n) models [U, as these are just a special case of our generalized 
multiple Ising model. 
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